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ABSTRACT 

The dynamics of Lennard- Jones fluid is studied through extended mode couphng theory 
(MCT) with the inclusion of the slow mode of defect density. Inclusion of defect density 
facilitates the liquid like state for temperatures much lower than predicted from ideal MCT. 
From the present model the temperature dependence of the power law exponent is obtained 
at a constant pressure. We have also computed the wave number dependence of the power 
law exponent. 
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Introduction 



The phenomena of Glass transition has been studied extensively through theoretical [|l], 
1^, ^ models, computer simulations H, 0] as well as scattering experiments p|, ||, |lOl- If a 
liquid is cooled beyond the freezing point avoiding the first order transition to an ordered 
crystalline state, it continues in the amorphous state with solid like properties. In this su- 
percooled state the liquid develops very long relaxation times and the transport coefficients 
(viscosity) abruptly increase. A microscopic mechanism for this slow dynamics using a sta- 
tistical mechanical approach has been proposed through the Mode coupling theory (MCT). 
This is based upon the idea that the transport coefficients in the supercooled state gets a 
feedback from the nonlinear coupling of slow modes in the fluid. The effect of the structure on 
the dynamics is studied from the wave vector dependence of the model. MCT predicts a two 
step relaxation process for glassy dynamics : power law (/5-relaxation) at intermediate time 
scales which crosses over to stretched exponential (a-relaxation) decay of density correlation 
function over long time scales. In the ideal model [l^ there is a sharp transition from er- 



godic to a nonergodic phase at a critical density. In the later version of more careful analysis 



1^, |T^, |15| the sharp transition is flnally removed from the ergodicity restoring mechanisms 
coming from the coupling of the density fluctuations with currents in the compressible fluid. 
The flnal ergodicity restoring mechanisms in a supercooled liquid behaves like a collective 
slow mode. There has been phenomenological extension of the theory to consider extra slow 
modes that develop in the liquid in the supercooled region. In the present work we consider 
an extension p!6| , 0, |18|, |T9[ of the simple feed back mechanism of ideal MCT through the 
inclusion of the extra slow mode of defect density in the hydrodynamic description for the 
isotropic liquid. Here we use this model to study the dynamics in the simple One component 
Lennard- Jones fluid in the supercooled state. In this respect we like to note that previous 
MCT studies on LJ liquid show that one component system freezes into a nonergodic 
phase at a relatively low density. 
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The defect density considered in the model can be associated with the vacancies in the 
sohds or void or free volumes in the fluid. The movement of these free volumes were consid- 
ered crucial for the transport in the liquid . Following Ref . |]l9l we consider the defects to 
be moving in a metastable potential well. Here the extent of coupling between defect density 
and the particle density governs the location of the ideal glass transition point that occurs 
as a result of the feedback mechanism proposed in MCT. Thus at a constant pressure we 
observe a transition line whereas a transition point is predicted in the simple MCT. Phase 
diagram of coupling parameters with temperature is studied through nonergodicity parame- 
ter. We compute the relaxation dynamics of fluid in the power law regime and calculate the 
power law exponents. From the present model the temperature dependence of the relaxation 
exponents arises naturally as a consequence of the critical line. The paper is organized as 
follows. In next section we briefly describe the model considered followed by section III 
where we discuss briefly the static strucutre factor used in the calculation. In section IV the 
results for the power law relaxation as predicted by the mode coupling theory is presented. 
In last section a brief discussion of the work is presented. 



In the mode coupling theory we focus on the normalized density-density correlation func- 
tion {'ipq{t)). The dynamical evolution of ipq{t) is written in terms of an integro-differential 
equation, 



where the constant = j3mS{q) / with m as the mass of the liquid particle and /3 is 
the Boltzman factor. S(q) is the static structure factor for the liquid. The bare transport 



The Model Studied 



%{t) + + / dsV^%t - s)i;q{s) = 0, 





(1) 



coefficient F^ |21| is defined in a dimensionless form as 




(2) 
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where q* = qa and jn{x) is the n order spherical Bessel function The Enskog time 



tE p3[ for the LJ fluid is approximated using the equivalent hard sphere diameter. The 



generalized transport coefficient Tq{t) has the bare (F^) as well as the mode coupling (r^'^(t)) 
contributions from the nonlinear coupling of density fluctuations. Following Ref. |jT9|, the 
mode coupling part of the generalized transport coefficient is given by, 

— * — * 

with ki = q — k. Here (pk-^it) is the normalized defect auto correlation function and gives the 
contribution to the transport coefficient coming from the coupling of the defect density with 
the density fluctuations. The mode coupling vertices with full wave vector dependence are 
given by 

k) = ^[2yU{q, k)S{k)S{ki) + Knoc{k)S{k)] (4) 



no 



and 



V^'\q, k) = ^[U'iq, k) - 2yU{q, k)]S{k)S{k^) (5) 
2no 



These equations are obtained |jT9| from the nonlinear equations for the extended set of 



hydrodynamic variables. The vertex function U in is given by, 

[/(g,A;) = ^[g.A;c(A;) + g.(A;i)c(A;i)] (6) 

where no is the particle number density. It is assumed that the defects are moving in a 
metastable double well potential with /t and y as two dimensionless parameters characterizing 
the coupling of the defect densities with the particle density and the depth of the potential 
well. Defect density is assumed to be weakly interacting with the mass density and is treated 
as a variable similar to mass density. For k = y = 0, (Bf) reduces to the standard result, ||TT|. 



The presence of the linear term in the mode coupling vertex comes through the coupling with 
the very slowly varying defect correlation In the next section we study the implications 
of this feed back mechanism that arises from the coupling of the density fluctuations as well 



as the defect density coupling. For this the static structure factor is taken as an input - 
measuring the extent of structural effects on relaxation. We use the standard model [^] for 
the Lennard- Jones Structure factor. We give a brief description of computation of the static 
structure factor next. 

Static Structure Factor 

In mode coupling equations interparticle interaction potential enters as an input via static 
structure factor S{q). For a Lennard- Jones (LJ) fluid the perturbation approach described 
in 1^ is being used for the computation of the structure factor. The full Lennard- Jones 
potential, 

U{r) = Mi^/rr - {^M'] (7) 
is divided into a repulsive part Uo{r) and an attractive part UA{r) as, 



f/(r) + e, r < Tq 

0, T > To 

J — e r < To 

I f/(r), T > ro 



with ro = A trial system with hard sphere diameter d < ro and an attractive potential 

energy part [/^(r) is taken. One then finds the appropriate diameter d through perturbation 
approach until e'^^^^^^ gxir) = e^'^^^^g{r). Where, subscript T represents the trial system. 
Once the hard sphere diameter d is known the static structure factor for the LJ liquid is 
approximated by. 



S{q) = STiq)+p I rfr[e^^-M(?T(r)][e- 



For the trial hard sphere system the Percus-Yevick solution with Verlet-Wiess correction 



271 is used. 



Power law exponents 

The ideal transition from an ergodic phase (fluid) to a non-ergodic (glass) phase is charac- 
terized in terms of non-ergodicity parameter defined as, tpq{t —>■ oo) = f{q)- In the long time 
limit or z — > limit of the Laplace-transformed form of the equation (|^) reduces to 



coupled equations are solved self-consistently for f{q). In solving the integral equations a 
maximum cutoff for the wave vector integration Aa = 30.0 with 200 grid points is used. A 
set of critical parameters (kc, i/c, Pc, T^) are identified below which all f{q)'s converge to zero 
values (ergodic phase) and above that critical point all f{q)'s converge to nonzero values 
(nonergodic phase). The critical points (t/c, ^^c^ Pc, in T — k phase space are located by 
fixing i/c = 0.15 and varying k, and T. Here we study the equations with the thermodynamic 
parameters changing in a way so that the volume change with the temperature keeping the 
Pressure constant, P* = Po^/t = 0.5 as shown in Fig. 1. Then, at a fixed pressure, we 
observe a line of transition giving different critical temperatures with different values of the 
metastability parameters. In the ideal model under similar conditions a single transition 
point is obtained. Along this transition line at a fixed pressure we calculate, from the 
solutions of the MCT equations, the power law exponent for relaxation in the LJ fluid as a 
function of temperature and density. We study the relaxation over time scales longer than 
the microscopic time scales t > (fi^r)"^, where r is the Lennard- Jones time (mcr^/e)^, but 
shorter compared to the a-relaxation time scales. The dynamical equation (0), for density 
correlation function ipq{t) is solved self consistently along the transition line shown in Fig. 
1. The metastable potential well for defects used for computing vertex functions is shown in 
Fig. 2 at critical parameters (Tc = 0.544, j/c = 0.15, = —0.038). The observed relaxation 





where vertices V"*-^-* and V^"^^ are given by equation (^ and equation (|]). These set of 
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behavior is given by, 

Mt) = fc{q) + Ch{q){t/to)-'' (10) 

where C is a constant and to is the time scale over which the power law relaxation persists. In 
Fig. 3 we show the decay of density correlation function with time (t/ r). The fit to the power 
law exponent is shown by dotted lines and the solid line represents the result from MOT 
equations for temperature = .544 in the power law regime. Temperature dependence of 
the power law exponent a is studied at a constant pressure P* = 0.5. The exponent a is 
calculated at the transition temperatures for fixed yc = 0.15 and the corresponding k'^s. Fig. 
4 shows the variation of the power law exponent a with temperature. We also calculate the 
wave vector dependence of the power law exponent a. Fig. 5 shows the exponent a along the 
wave vector go" at T = 0.544 and p = 0.958. In the inset the static structure factor is shown 
to depict the alignment of the maxima of S{qa) with minima of exponent a and vice- versa. 
The arrows are shown at the maxima and minima. The power law exponent is minimum 
at the peak of the static structure factor, which reflects the characteristic slow relaxation at 
the diffraction maximum for liquids. 

Discussion 

In the simple mode coupling theory there is a single transition point with a fixed T^. As a 
result of which the power law exponent determined in its vicinity is temperature independent. 
In the present case through the coupling of the density fluctuations with defect motions a 
line of dynamic instability is predicted and this results in a temperature dependent exponent. 
Here we consider the simple one component LJ system and we obtain an oscillatory behavior 
of the power law exponent a - as shown in Fig. 4. For hard sphere system a similar oscillatory 



behavior is found [28| from the solutions of the MCT equations. This was also shown to 



agree with the power law exponent that was calculated from fitting of the light scattering 



data on Colloids |T0|. A similar behavior was observed in the wave vector dependence of the 



exponent as well. Computer simulations of one component systems at metastable densities 
are particualrly rare due to problem of crystallization and at this time we are unable to find 
simulation data to look for such behavior. However simulation results for binary mixtures 



exist and we studied the data reported in ref. without restricting to a fixed non-ergodicity 
parameter (NEP) as in the ideal model that gives rise to a single exponent. We do observe 
in this case a qualitatively same oscillatory behavior for the exponent a of the simulation 
data as shown in Figure 6. The NEP that is obtained in our fitting also shows qualitatively 
similar behavior as predicted from the extended model studied here. The present model 
can be extended to the a-relaxation regime as well where the combined role of the density 
fluctuation and the defect fluctuation need to be considered. 
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Figure Captions 

Fig 1 : Phase diagram in T—k space at = 0.15 along a constant pressure line. Temperature 
is expressed in units of e/kB- 

Fig 2 : Variation of depth of the potential with n/n* for T* = 0.544. In the figure h*[n) 
represent the dimensionless quantity h{n)f3ena^. 

Fig 3 : Decay of normalized density-density correlation function ipglt/r) with time (t/r) at 
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T = 0.544. Where r is the Lennard- Jones time (mcT^/e)^. Dotted hnes show the power law 
fit (|10D (see text). 

Fig 4 Variation of power law exponent a as a function of temperature T. T is expressed in 
units of e/fcs. 

Fig 5 : Variation of exponent a as a function of wave number qa for a packing fraction 
T = 0.544. In inset static structure factor is shown along wave number qa. 



Fig 6 : Variation of power law exponent a with temperature as obtained by fitting [lO|(see 



text) to the self correlation function data from Ref. p9[ at wave vector q = 7.22. 
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